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I. INTRODUCTION 

There was recently much activity in the phenomenol- 
ogy of hadronic heavy quark pair production in connec- 
tion with the Tevatron and the CERN Large Hadron Col- 
lider (LHC) which will have its startup this year. There 
will be much experimental effort dedicated to the discov- 
ery of the Higgs boson. There will also be studies of the 
copious production of top quarks and other heavy par- 
ticles, which also serve as a background to Higgs boson 
searches as well as to possible new physics beyond the 
standard model. Therefore, it is mandatory to reduce 
the theoretical uncertainty in phenomenological calcula- 
tions of heavy quark production processes as much as 
possible. 

Several years ago the next-to-next-to-leading order 
(NNLO) contributions to hadron production were cal- 
culated by several groups in massless QCD (see e.g. [l[ 
and references therein) . The completion of a similar pro- 
gram for processes that involve massive quarks requires 
much more dedication and work since the inclusion of an 
additional mass scale dramatically complicates the whole 
calculation. 

Until very recently there was the belief that the next- 
to-leading order (NLO) description of heavy charm and 
bottom production in hadronic collisions considerably 
underestimates the experimental results. In recent, more 
refined analyses 0, SHI it was shown that a NLO anal- 
ysis does in fact properly describe the latest charm and 
bottom quark production data ||. The authors of @ 
and [1, 0] deal very differently with the problem of large 
mass logarithms which constitute the central problem in 
the heavy quark phenomenology. Data on top quark pair 
production also agrees with the NLO prediction within 
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theoretical and experimental errors (see e.g [6(). In all 
of these NLO calculations there remains, among others, 
the problem that the renormalization and factorization 
scale dependence of the NLO calculations render the the- 
oretical results quite uncertain. This calls for a NNLO 
calculation of heavy quark production in hadronic colli- 
sions which is expected to considerably reduce the scale 
dependence of the theoretical prediction. 

At the lower energies of Tevatron II, top quark pair 
production is dominated by qq annihilation (85 %). The 
remaining 15% come from gluon fusion. At the higher 
energy LHC, gluon fusion dominates the production pro- 
cess (90 %) with 10 % left for qq annihilation (percentage 
figures from [6j). This shows that both qq annihilation 
and gluon fusion have to be accounted for in the calcula- 
tion of top quark pair production. 

In general, there are four classes of contributions that 
need to be calculated for the NNLO corrections to the 
hadronic production of heavy quark pairs. The first class 
involves the pure two-loop contribution, which has to be 
folded with the leading order (LO) Born term. The sec- 
ond class of diagrams consists of the so-called one-loop 
squared contributions (also called loop-by-loop contribu- 
tions) arising from the product of one-loop virtual matrix 
elements. This is the topic of the present paper. Further, 
there are the one-loop gluon emission contributions that 
are folded with the one-gluon emission graphs. Finally, 
there are the squared two- gluon emission contributions 
that are purely of tree-type. 

Bits and pieces of the NNLO calculation are now be- 
ing assembled. The recent two-loop calculation of the 
heavy quark vertex form factor [3] can be used as one of 
the many building blocks in the first class of processes. 
In this context we would also like to mention the recent 
work [8j on the NNLO calculation of two- loop virtual am- 
plitudes performed in the domain of high energy asymp- 
totics, where the heavy quark mass is small compared 
to the other large scales. In this calculation mass power 
corrections are left out, and only large mass logarithms 
and finite terms associated with them are retained. The 



2 



authors of the present paper have been involved in a sys- 
tematic effort to calculate all the contributions from the 
second class of processes, the one-loop squared contribu- 
tions. We shall describe the present status of this pro- 
gram in the next paragraph. In the work [9( the full, 
exact NLO corrections to tt+jet are presented. When in- 
tegrating over the full phase space of the jet (or gluon), 
this calculation can be turned into a NNLO calculation 
of heavy hadron production of the third class. To our 
knowledge there does not exist a complete calculation 
of the fourth class of processes, the squared two-gluon 
emission contributions. 

Let us briefly describe the status of our effort to cal- 
culate the one-loop squared contributions for the second 
class of processes. The highest singularity in the one-loop 
amplitudes from infrared (IR) and mass singularities (M) 
is, in general, proportional to (1/e 2 ). This in turn implies 
that the Laurent series expansion of the one-loop ampli- 
tudes has to be taken up to 0(e 2 ) when calculating the 
one-loop squared contributions. In fact, it is the 0{e 2 ) 
terms in the Laurent series expansion that really com- 
plicate things [nj since the 0(e 2 ) contributions in the 
one-loop amplitudes involve a multitude of multiple poly- 
logarithms of maximal weight and depth 4 [ll| . All scalar 
master inte gral s needed in this calculation have been as- 
sembled in [id [ll| . Reference [1(3] gives the results in 
terms of so-called L functions, which can be written as 
one-dimensional integral representations involving prod- 
ucts of log and dilog functions, while [ll[ gives the results 
in terms of multiple poly logarithms. The divergent and 
finite terms of the one-loop amplitude qq — > QQ were 
given in The remaining O(e) and 0(e 2 ) amplitudes 
have been written down in Squaring the one-loop 

amplitudes leads to the results of the present paper. In a 
recent work fl4| we have presented closed-form, one-loop 
squared results for heavy quark production in the fusion 
of real photons. 

In this paper we report on a calculation of the NNLO 
one-loop squared matrix elements for the process qq — > 
QQ. The calculation is carried out in the dimensional 
regularization scheme [15| with space-time dimension 
n = 4 — 2e. In sequels to this paper we shall present 
results on the square of hadroproduction amplitudes orig- 
inating from the gluon fusion subprocess gg ~ > QQ and 
photoproduction amplitudes 73 — ► QQ. 

In our presentation we shall make use of our notation 
for the coefficient functions of the relevant scalar one- 
loop master integrals calculated up to 0(e 2 ) in [l(|. For 
the case of gluon-gluon and quark- antiquark collisions, 
one needs all the scalar integrals derived in e.g. the 
one scalar one-point function A; the five scalar two-point 
functions B\, B2, -B3, B4, and B$; the six scalar three- 
point functions C%, C2, C3, C4, C5, and Cq; and the three 
scalar four-point functions Di,Z?2, and D3. Taking the 
complex scalar four-point function D2 as an example, we 
define successive coefficient functions for the Lau- 



rent series expansion of £>2- One has 

D 2 = l C e {m 2 ){-A- 2) + Id™ + + eD^ 

+e 2 D^+0(e 3 )}, (1.1) 

where C £ (m 2 ) is defined by 

We use this notation for both the real and the imaginary 
parts of Z?2, i-e. for RcZ?2 and ImZ^. Similar expansions 
hold for the scalar one-point function A, the scalar two- 
point functions Bi, the scalar three-point functions Ci, 
and the remaining four-point functions Di. The coeffi- 
cient functions of the various Laurent series expansions 
were given in jlG ] in the form of so-called L functions, 
and in [llT ] in terms of multiple polylogarithms of max- 
imal weight and depth 4. It is then a matter of choice 
which of the two representations are used for the nu- 
merical evaluation. The numerical evaluation of the L 
functions in terms of their one-dimensional integral rep- 
resentations is quite straightforward using conventional 
integration routines, while there exists a very efficient al- 
gorithm to numerically evaluate multiple polylogarithms 

Let us summarize the main features of the 
scalar master integrals. The master integrals 
A, B±, B3, -B4, C2, C3, and -D3 are purely real, whereas 
B 2 , B 5 , Ci, C4, C5, Cq, Di, and D 2 are truly complex. 
From the form (AB* + BA*) = 2 (Re A ReB + ImA ImB) 
it is clear that the imaginary parts of the master integrals 
must be taken into account in the one-loop squared con- 
tribution. The master integrals B2, -B5, Ci, C4, C5, and 
Cq are (t ^> u) symmetric, where the kinematic variables 
t and u are defined in Sec. [TTJ 

The paper is organized as follows. Section |TT] con- 
tains an outline of our general approach and discusses 
renormalization procedures. Section IIIII presents NLO 
results for the quark-antiquark annihilation subprocess. 
In Sec. IIVI one finds a discussion of the singularity struc- 
ture of the NNLO squared matrix element for the quark- 
antiquark annihilation subprocess. In Sec. [V]we discuss 
the structure of the finite part of our result. Our results 
are summarized in Sec. IVII In the Appendix we write 
down expressions for the building blocks of that part of 
the finite result that originates from the square of box 
diagrams. 

II. NOTATION 

The heavy flavor hadroproduction proceeds through 
two partonic subprocesses: gluon fusion and light quark- 
antiquark annihilation. The first subprocess is the most 
challenging one in QCD from a technical point of view. It 
has three production topologies already at the Born level. 
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FIG. 1: The lowest order Feynman diagram representing light quark-antiquark annihilation. Normal solid lines represent the 
light quarks, the curly line represents the gluons and the thick solid lines correspond to the heavy quarks. 



Here we consider the second subprocess where there is 
only one topology at the Born term level (see Fig.[l}. Ir- 
respective of the partons involved, the general kinematics 
is, of course, the same in both processes. In particular, 
for the quark-antiquark annihilation, Fig. [TJ we have 

q(pi) + q{p2) ^Q(p 3 ) + Q(p4), (2.1) 

The momenta directions correspond to the physical 
configuration; e.g. p\ and p 2 are ingoing whereas p 3 and 
P4 are outgoing. With m being the heavy quark mass, 
we define 

s = (pi +p 2 ) 2 : t = T-m 2 = (pi-p 3 ) 2 -m 2 , 

u = U-m 2 = (p 2 -p 3 ) 2 -m 2 , (2.2) 

so that the energy-momentum conservation reads s + 1 + 
u = 0. 

We also introduce the overall factor 

C = (giC e {m 2 ))\ (2.3) 

where g s is the renormalized strong coupling constant 
and C e {m 2 ) is defined in \l.2\ . 

As shown e.g. in (l2l. 1 1 3t| the self-energy and ver- 
tex diagrams contain ultraviolet (UV) and infrared and 
collinear (IR/M) poles even after heavy mass renormal- 
ization. The UV poles need to be regularized. 

Our renormalization procedure is carried out as fol- 
lows: when dealing with massless quarks we work in the 
MS scheme, while heavy quarks are renormalized in the 
on-shell scheme, where the heavy quark mass is the pole 
mass. For completeness we list the set of one-loop renor- 
malization constants that we have used: 

Z i = 1 + ~\ H N c - ni)C £ (n 2 ) - C £ (m 2 )} , 
Z m = l-g 2 C F C E (m 2 ) £ 3 { -_ 2 ^ y 

Z 2 = Z m , (2.4) 
a 2 

Z 1F = Z 2 ~^N C C E (^ 2 ), 

£ 

Z lf = l-^N c C e (ti 2 ), 

Z 3 = 1 + 9 f {( 5 -N c - \ni)C s {p?) - \c e {m 2 )) 



= 1 + 9 ~f {(#> - 2N c )C e (p 2 ) ~ \c e {m 2 )} , 

Z « = 1 -T {y^(y)-^ £ (m 2 )}. 

with 0o = (UNc — 2ni)/3. ni is the number of light 
quarks, Cf = 4/3, and Nc — 3 is the number of col- 
ors. The arbitrary mass scale \x is the scale at which 
the renormalization is carried out. The above renormal- 
ization constants are as follows: Z\ for the three-gluon 
vertex, Z m for the heavy quark mass, Z 2 for the heavy 
quark wave function, Z\p for the (QQg) vertex, Z\f for 
the (qqg) vertex, Z$ for the gluon wave function, and Z g 
for the strong coupling constant a s . Note that Z\ is not 
actually needed in the present application, but we have 
presented it for completeness. For the massless quarks 
there is no mass and wave function renormalization. 

The above coefficients (except for Z g ) are needed if 
one renormalizes graph by graph. However, one could 
choose another route. From the field-theoretical point of 
view, the renormalized matrix element is obtained from 
the unrenormalized one by 

M Ien — Y\ Z c I -Mbare(flWe — * Z g g s , TObaxe Z m m r ), 
n 

where Zt n are the wave function renormalization con- 
stants for all the external on-shell particles under con- 
sideration. If one formally expands Mbare (e.g. -Mbarc = 
^-^o + g 2 ^h + • • ■) and the renormalization parameters 
Zf n as a series of powers in the coupling constant to the 
requisite order, one arrives at the one-loop order result 

^ ^ 1/2 

Ml <Ten = ]J n Zj M ((7baro — * Z g g s , mbare —* Z m m r ) 
+g 2 M 1 {g s ,m r ), (2.6) 

where now the Zt nt \ correspond to the one-loop renor- 
malization constants for the external particles. In our 
case one has Zf lt i = Zt 2t i = 1 and ^/ 3 ,i = Z/ 4 ,i = Z 2 - 
Thus, one could apply inverse wave function renormaliza- 
tion for external legs and then replace the bare coupling 
constant <?barc — ► Z g g s (as the mass parameter m does 
not explicitly enter the leading order Born term matrix 
element, it is not renormalized at that order). We have 
verified that, in both ways, we arrive at the same renor- 
malized result. 
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FIG. 2: One-loop Feynman diagrams contributing to the subprocess qq — > QQ. The loop with the dotted line in (g) represents 
the ghion, ghost, and light and heavy quarks. 



In order to fix our normalization we write down the 
differential cross section for qq — > QQ in terms of the 
squared amplitudes \M\ 2 . One has 

da ^ ^ = ^kJ_|M| 2 _ -, (2.7) 
2s 47V,?, 1 'w-W v ; 



qq^QQ 



where the n-dimensional two-body phase space is given 
by 



d(PS) 2 



-2s 



(47T) 



8?rs r(l-e) 



5(s + t + u)dtdu 
. (2 ' 8) 

and we explicitly show flux (Ap\p2) 1 = (2s) , initial 
quark and antiquark spin (2s f + I) -2 = 1/4, and color 
Ng 2 averaging factors. Then, at the leading Born term 
order for qq — > QQ, we have 

1 
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4 ,^Ilo = 16 



t 2 + u 2 



B. 



(2.9) 



III. NEXT-TO-LEADING ORDER RESULT 



Folding the one-loop matrix elements depicted in 
Fig. [5] with the LO Born term, Fig. [TJ one obtains the 
virtual part of the NLO result. Although NLO virtual 
corrections to heavy flavor hadroproduction were calcu- 
lated before for the qq — > QQ case, one cannot find ex- 
plicit results for this subprocess in the literature. We 
have therefore recalculated the virtual NLO contribution 
to gg-annihilation. In fact, we have calculated the virtual 
NLO results up to 0(e 2 ). As it turns out, the expressions 
for the NLO virtual e 1 and e 2 contributions considerably 
simplify the presentation of the corresponding NNLO re- 
sults, in as much as they appear as important building 
blocks in the NNLO results. 
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First, we write down a few abbreviations that we shall use throughout the paper: 



= y/1 — 4m 2 /s, D = m 2 s — tu, 
z% = s + 2t, Z2u = s + 2u, z t — 2m 2 + t, z u — 2m 2 + u. 

Note that D is not the space-time dimension. We further define the functions: 
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■\C F - ^N C + 1 - ^ 



B{ 3) I 3(7 



5N C 2m N c 
3 3 2/3 2 



) + C^Nc 7 ^ - {c[ 3) s ~ cf (2m 2 - a)} (2C F - N c ), 



(3.1) 



(3.2) 



F^> = 2 (s/3 2 (2C F ~ N c ) - l2m 2 N c ) - B 2 3) s0 2 (2C F - N c ) + B ( 5 3, (8m 2 + s)N c + C\ 3} Qm 2 sN c , (3.3) 



= — {2 



5G 



+e 



1 



s„r ( - - 

CO 



%2 



t s 2 [3 2 
At 2 



s/3 2 



s s 2 /3 2 ! ' • 



-Ci 3) 2t [1 + 2- 



T 



-8m' 



3 _ 2z2_ 

t s 2 /3 2 



BY '2 



B i3) 2 12 



^i(i) I r, 2 , 8m 2 t + s 
-Cf J 8m 2 + 4s H 



s/3 2 



J 2 m 2 s 2 + 2t 3 



(CP-Dg't) [ 3s + 4t + 2st 



,C0 



s - i 



+e z 



,3^ 



.0) 



2/ 2 



^COv 



C0 + 2 



(c| j) -^^)^(2^ + s 2 + 2t 2 ) 



>C0v 



g2 
S^ 2 



Ci i} 2- [ s-4i-2st- 



t 



s - t 



D 



(3.4) 



F 



0) 



8m 2 [ - 

u 



Z2u 

s 2 (i 2 



D 



B 



CO 2z t 
sf3 2 



8m 4 



S ~ ~ - C«2zt ( 1 + 2 U - ) + (cf - D 2 W«)i(2i3 + s 2 + 2zt 2 ) 



u 

CO 



2^2^/ 

I 2 /? 1 



s 2 /3 2 

It 



(3.5) 



2u 
s 



v\:m---)-B^2[i-^ 



fco , 



4z„ + 



8m 2 u 



s/3 2 



9 2 

m s 



c^^ + ci^ + c^-^v 



2m 2 



D 

•<C0„ 



o 2 Z 2u 

2m s-^j- 



4u 



2z, 



4u 



stu 



The additional subscript "u" in some of the scalar coef- 
ficient functions in the expression for F^ (such as B[ 3 J) 
is to be understood as an operational definition prescrib- 
ing a (t u) interchange in the argument of that func- 
tion, i.e. = B^l^, etc. Note that B 3 5 , C{, and 
C\ are intrinsically (t <-* u) symmetric (see [13]). Taking 
the (t <-> u) symmetry of B§, C{, and C\ into account, 
one notes a corresponding (t <-> zt) symmetry for the first 



and third square brackets in F^ and F 4 . 

Before presenting our result for the NLO matrix el- 
ement, we would like to comment on its color struc- 
ture. First note that all the vertex and self-energy (VSE) 
graphs are proportional to the LO Born term color ma- 
trices (see Refs. [IUdlj])- Both the parallel ladder box, 
Fig. [211a), and the crossed ladder box, Fig. (Hb), have 
their own color structures. Altogether one obtains the 
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following three color structures, 



in W$" represent the contributions from boxes a 
and b, respectively. 

The massless limit of our NLO result Eq. (|3.7j) without 
the 0(e) and 0(e 2 ) order terms was compared (includ- 
ing also the imaginary part) with corresponding results 
obtained from the methods developed in Ref. [TtJ [22| - 
There was agreement . This serves as a rigorous check 
on our singularity structure as well as on all the mass log- 
arithms of our original NLO matrix element [l2l |. 



r<Sl 



Tr(T a T b ) Tr(T a T b ) 
r Fi(T a T* } T c } Tr(T"^T' a T lC ) 



N c 



2 



(3.6) 

7 

3' 



' 4 N c 



from folding the Born term with the VSE graphs, the 
parallel ladder box, Fig.[^a), and the crossed ladder box, 
Fig. [D(b) , in that order. The common factor o!a = N c — 
1 = 8 is the dimension of the adjoint representation of 
the color group SU (Nc)- We present our NLO result 
separately for these three color structures. 

At NLO the final spin and color summed matrix ele- 
ment can be written down as a sum of three terms: 



IV. SINGULARITY STRUCTURE OF THE 
NNLO SQUARED AMPLITUDE 

The NNLO final spin and color summed squared ma- 
trix element can be written down as a sum of five terms: 



1 2 
~5 Try I ^ I Loop X Born 



Re 



\w^(s) + -W^ 1 \e) 



+W^(e)], 



(3.7) 



where C has been defined in (|2.3p . The notation 
I M | Loop x Bom means that one is retaining only the O(o%) 
part of |Af| 2 . 

The first two coefficient functions in (|3.7[) have a rather 
simple structure: 



W { - 2) {e) = -25(2C F - N c + 3) 
W (_1) (e) = -2B(5C F 

x (2C F 



(3„ 



°4 s ~ °6 



(2m 2 - s) 



9 r — t —n 
Nc)- - 71n(— ) + 21n(— ) 
3 L m m 



where B is the Born term defined in Eq. (|2.9|) . One 
should keep in mind that the overall Born term factor B 
above contains a term multiplied by e. Therefore, if the 
expression for B, Eq. ||235D, is substituted in W^ 2) and 
W^~ x ' , we will obtain (e) -1 and finite terms from the 
first two terms of Eq. (|3.7[) . 

The third term in Eq. (|3~7|) reads 



NLO' 



where 



and where 



NLO 



(i) 



W, 



U) 



w. 



0) 



(3.9) 



(3.10) 



Re 



ly(-4) (£)+ ly(-3) (£) (41) 



+ ±V^ 2 He) + -V^ 1 \e) + V^(e) 



q I M | Loop x Loop 



where C has been defined in (|2.3p . Note Eq. (|4.ip is not 
a Laurent series expansion in e since the coefficient func- 
tions V {m) (e) are functions of e as explicitly annotated 
in Eq. (|4.1[) . It is nevertheless useful to write the NNLO 
one-loop squared result in the form of Eq. (|4.1|) in order 
to exhibit the explicit e structures. All five coefficient 
functions V {m \e) are bilinear forms in the coefficient 
functions that define the Laurent series expansion of the 
scalar master integrals Some of these coefficient 

functions are zero and some of them are just numbers or 
simple logarithms. In the latter case we will substitute 
these numbers or logarithms for the coefficient functions 
y( m ) j n the five terms above. This will be done for the 
coefficient functions A^ m \ B { ~ 1] , b{~ 1] , 



(-i) 

3u 



and C, 

We mention that in the course of our work we took full 
advantage of the fact that in [l2| all the poles in the ma- 
trix element for the qq — > QQ subprocess are multiplied 
only by the leading order Born Dirac structure to cast 
the singular terms of the squared matrix element into an 
appropriately factorized form. 

Before proceeding further, we present three more color 
structures appearing in the NNLO calculation in addition 
to the ones presented in Eq. (|3.6[) : 



W[ i] = 2B F[ j) 



771 _ 

128— —Fo 



U) 



W ( 2 3) =-2B/? ln 1+J ' T " 



s 4 /? 4 



W. 



F 



U) 



(3.11) 



Note that the first term in (|3.11|) originates entirely from 
the sum of self-energy and vertex diagrams while the sec- 
ond term is due to renormalization. The terms FP' and 



brparjia ry\b \ 



16 



Ni 



N 2 C 



19 



i 



N 2 



d,A 
16 



N 2 



(4.2) 
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The above three color structures arise from folding box a 
with box a, box b with box b, as well as the interference 
of the two boxes, respectively. 

Let us first introduce a notation which will help us to 
present the coefficients of the singular terms in the most 
concise fashion: 



L x = (2C F - N c ) [Ct 1] s ~ Cf x >(2 

—t —11 
L 2 = 15C F - 141n( — ) - 41n(— ), 
in m z 
—t — 11 

L 3 = 35CV - 381n(— -41n(— =■), 

— t —u 
U = 5C F - 21n(— ) - 41n( — ). 



•»(-!)/ 



«))> 



(4.3) 



The two most singular terms in (|4.1[) are proportional 
to the Born term B defined in (12.91). One has 



V { ~ 4 \e) = (2C F - N c + 3) 2 B, 
V ( - 3 \e) = 2{2C F - N c + 3)B 

We also obtain 



(4.4) 



V^ 2 \e 



B r 

~3 



(3L 1 +L 2 )(L 1 +5C f )* (4.5) 



21n(— )(7Li + L 3 ) 



-41n(^)(L 1+ L 4 ) 



-{2C F -N c + Z)F l 



(o) 

NLO ' 



The last term in Eq. (|4.5p is obtained from folding the 
0(e~ 2 ) singular term of the matrix element with its finite 
part, while the rest is obtained from folding the single 
poles. Note that when one substitutes the Laurent ex- 
pansions for B and F^ , one gets additional 1/e poles 
and finite terms in Eq. (|4.5|) . 



The structure of the last term in Eq. (|4.1|) is a little 
more complicated. One has 



v<-'> (e ) = -l;< 



(0) 



w. 



(oh 



(4.6) 



-—Fl 
7 a 

-(2C F 



(o) 



(0) 



N c + 3) 



NLO 



V 



The terms multiplied by the L m functions above are due 
to folding the single pole terms in the original matrix 
element with its finite O(e ) part, while the last term is 
due to interference 0(e~ 2 ) x 0(e) terms in the original 
matrix element. This pole term is due to the Taylor 
expansion of the original matrix element and cannot be 
deduced from the knowledge of the LO terms alone. The 
function i^iio ^ s defined in Eq. (|3 . 1 Ojl and is nothing 
but the finite part of the NLO term with indices of the 
coefficient functions of the scalar master integrals and the 
power of the logarithm that multiplies the (3q function, 
shifted upward by 1. For the remaining term V , one 
obtains 



V = -2B 



A, 2 y, i(! „ N c 2rn + 2 + 28N c , d( o) 2(3 2 - lggg + 9N C 
_l n( _) +8 ^-_ ^ + B, 



27 



9 



(4.7) 



+^ 0) ^(5iV c + n / -9C F ) 



56 

Y 



1 



Z2 
2/92 



8m' - 



t s 2 f3 
2z 2 



t s 2 (3 2 

2 z 2u 



m 2 D 
2 s-t 

s + — 

8 7s - 41 

s 



2(<o0 2 C F - N c ) - B^ ] 2(3 2 {2C F - Nc) - B^''2N C - C^'sNc 



(0), 



D 



[ci 0) S z t + Ci 0) 2t 2 + C ( fst-D^st 2 ) 



C[ a) sz t + C^2t 2 + Ci 0) s t-D^ s t 2 



s 2 (3 2 



-(-1) »^ c{° (-('• 



su 



SZ 2 i 

D 



2zt_ 

(3 2 



,(0) 

3m 

1_ . 2z2m 

u s 2 (3 2 



(0), n (o) 

s ~ u 2u 

8 9s + Au 



D 



' sz u 



-C^2u 2 + C^ S u-D^su 2 



When one substitutes the Laurent expansions for F^ , Laurent series expansion up to the finite term, these O(e) 
F 4 (0) , and F^Q, one gets finite and O(e) order terms in contributions should be omitted. 
Eq. (|4.6p . However, since we are only interested in the 
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V. STRUCTURE OF THE FINITE PART 

In this section we present the finite part of our result. 
We calculate the finite part in several pieces, e.g. 



W°> = Re 



V, 



(0) 



Bfi 



V, 



(0) 



Bh 



(0) 



/o/o 



(5.1) 



The first two terms originate from the interference of 
the 0(£ _1 ) x 0{e) and 0(e" 2 ) x 0{e 2 ) pieces of the ini- 
tial matrix element. Each of them can be conveniently 
presented as a sum of five compact expressions: 



v bI - G x + G 2 + G 3 + G 4 + G 5 , 



(5.2) 



where 

Gi = -128m 2 D(Ll + L 2 /3) 



F, 



(i) 



+ l2s(3 2 C F - 2siV c - B ( 2 ) 2s0 2 {2C F - N c ) 



- B { °'2sN c - C^'s'N c \ /(s 4 /? 4 ) , 

G 2 = -2B{L\ + L 2 /3) hlF^ - 2n t - 2 

- 28N C + 216G F - 27N C /P 2 - B 2 ;o) 3(18G F 



(°) „2 ; 



9N C - 2(i 2 ) - B { 5 0) 6(9C F - 5N C - m) /27 



G 3 = ^B\n 2 {—){L 1 +L 2 /3) 



It' 



(5.3) 



G, 



-16(7i* + L 3 ) F 3 (1) 3/112 
+ 8m 2 (lA - z 2 /(s 2 /3 2 )) + (Cf ^ + G^ 0) 2t 2 /s 
+ G| 0) t - D ( 2 k 2 ){2D + s 2 - st)/D~\ /3 , 

-32(LJ + L 4 ) [F 4 (1) 3/32 

+ 8to 2 z 2u /( S 2 /3 2 ) + B<°>2(2.D/(«0 - 1) 

- Bf ] 2z 2u /{sP 2 ) - c[ 0) (m 2 sz 2u /D 

- 2(8m 4 + st)/{s(3 2 )) - {C^2u/s + cf ] 
-D^u)(z 2u -stu/D)]/3. 



G s = 



The first three terms above are due to the VSE contri- 
butions, and the last two terms are due to the two box 
diagrams. Similarly, for the second term in Eq. (|5.ip we 
write 



with 



Hi = 
H 2 = 

H 3 = 
Hi = 

H* = 



V b% =H l + H 2 + H 3 + H i + H 5 , 



128(2C F -N c + 3)Dm 2 F 2 {2) + 4s/? 2 (7G F + Nc) - 10sN c - B i 2 1) 2s(3 2 (2C F - N c ) 



B^2sN c ~C^ S 2 N c }/{s i (3 i ) , 



~{2C F -N c + 3)5 F 1 (2) 162 + 2(1296G F + 76N C - 10nj - 10 - 2A3N c /0 2 ) + B y 2 >20 2 
- ^ 1} 18(18G F - 9iV c - 2/3 2 ) + Bf ] l2{N c + 2n t ) - B^ } 36{9C F - 5N C - n/)j/81 , 

2 

(2G F - N c + 3)B/3 ln 3 (^)/3 , 

-112(2C F - N c + 3) \fP 3/112 + 24to 2 (1/£ - z 2 /(s 2 /3 2 )) + (z t (2C[ 0) + c[ 1] ) 
+ 2t 2 {2Cf ] + Cf ] )/s + t{2Cf ] + cf } ) - t 2 (2D ( 2 0) + D ( 2 ] )){2D + s 2 - st)/D~\ /3 , 

-32(2G F - N c + 3) [F 4 (2) 3/32 + 8m 2 (l/u + z 2u /(s 2 f3 2 )) + B^2(2D/{su) - 1) 



>(0)r 



B<i ] 2z 2u /(s0 2 ) - (C?>z u + C^2u 2 /s + G 4 ' 



.(o), 



(°)o„,2, 



•<(o)„ 



D { 2 Ju 2 ){AD + 3s 2 + 2su)/D 



C[ 1) (m 2 sz 2u /D - 2(8m 4 + st)/(s/3 2 )) - (C$2u/a + g| X) - D ( 2 Ju)(z 2u D - stu)/D] /3 



(5.4) 



(5.5) 



Note again that the 0(e) and 0(e 2 ) order terms in the above expressions for Vgl and Vgl can be disregarded. 

We also mention that the scalar coefficient functions with 
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superscript "2" above involve multiple poly logarithms. 

We emphasize that the factorized forms of all the ex- 
pressions given in this paper hold only when one retains 
the full e dependence in the Born and NLO terms. 

The last term in Eq. (|5.ip comes from the square of the 
O(e ) term of the matrix element. It can also be written 
as a sum of five terms: 

V Mo = MvSE + m bvse + M aa + M ba + M bb . (5.6) 

The first term is the square of the finite parts of vertex 
and self-energy graphs; the second one is an interference 
of the vertex and self-energy graphs with the two box 
diagrams. These two terms can be presented in a very 
compact form: 



Ah 



VSE 



F, 



(o) 



W. 



(0) 



BR 



|F 2 (0) | 2 32m 2 J D/(s 5 /3 6 



(o) 



(5.7) 



Tf) 

-A)ln( — ) 
Mbvse = 7P + 2P| tM , 



W, 



f/2- 
F (o) 



w[ 0) -2BF[ 0) 



x (F 3 W +F 4 (0) ) 



(5.8) 



with 



P = 64m 2 F 2 (0) * 



2D/t - B[ 0) D/t 



c[ 0) Tz 2 ~ cf ] 2tT 



(C| 0) -D%"t)(D + ?) /(3s a /3 4 ) 



l(0) 



(5.9) 



When writing out P\t<->u one has to do the t *-> u op- 
eration in all the terms in the function P, i.e. for 

z,t,F<°\B[°\(f\T, andDf (cf and Cf } are t <- u 
symmetric) . 

Finally, we are left with the last three terms in 
Eq. (|5.6[) , which are the longest terms in our NNLO re- 
sult. However, to our surprise, we were able to discover 
nice factorization properties of the square of the two box 
diagrams. This part of the cross section can be put to- 
gether with the help of several building blocks; e.g. each 
of the last three terms in Eq. (|5.6[) can be written as a 
sum of bilinear products. Each of the factors in the bi- 
linear products are linear combinations of scalar integral 
coefficient functions multiplied by some combination of 
kinematic variables. To be more specific, we write 



76 

Maa = — 



+ m 2 Q 5 Q* 11 - 2s- 1 Q 6 Qt 2 + Q 7 Q 
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M ba 



-^g M a(1 \ t^U 1 

16 r 

— s~ 1 QsQ* 14 + 4:m 2 Q 9 Q* 15 + Q W Q* 16 



(5.10) 



+ 2m 2 Q 11 Q* w + 2s- 1 Q 12 Q* l7 + Qi 3 Q* ls 



Explicit expressions for the 18 linear forms Qi are given 
in the Appendix. The bilinear forms above arise from 
folding certain pairs of Dirac structures in our original 
matrix element. The expression for M ba represents the 
result of the interference of the finite parts of box a and 
box b. 

It is quite obvious that the factorized forms for the 
finite part of the NNLO result for the qq — > QQ sub- 
process should also hold for the corresponding massless 
amplitudes. We have not seen this being displayed in the 
literature. 

In the finite contribution, Eq. (|5.ip . one can see the 
interplay of the product of powers of e resulting from the 
Laurent series expansion of the scalar integrals [cf. Equa- 
tion (II. ip ] on the one hand, and powers of e resulting from 
doing the spin algebra in dimensional regularization on 

the other hand. For example, for the finite part one has 

( — 1) (o)* 

a contribution from Cg B\ ' as well as a contribution 

from C^B^*. Terms of the type C^ 1} b[ 0) *, where 
the superscripts corresponding to e powers do not com- 
pensate, would be absent in the regularization schemes 
where traces are effectively taken in four dimensions, i.e. 
in the so-called four-dimensional schemes or in dimen- 
sional reduction. 

We want to emphasize that all our factorized results 
given in this paper take up about 10 Kb of hard disk 
space. This has to be compared with the length of the 
original computer output. The original computer output 
for the one-loop squared cross section of the qq — > QQ 
subprocess turned out to be very long and took up ap- 
proximately 4 MB of hard disk space. Therefore, the 
reduction is of the order of 400 in the present case. 

As a final remark we want to emphasize that we have 
done two independent calculations using REDUCE [2(| 
and FORM [21( when squaring the one-loop amplitudes. 
After casting the results into the above compact form, 
we have checked the final result against the original un- 
treated versions using again the REDUCE Computer Al- 
gebra System. 



VI. 



CONCLUSIONS 



We have presented analytical 0(a A s ) NNLO results for 
the one-loop squared contributions to heavy quark pair 
production in quark-antiquark annihilation. These are 
the first exact results for the hadroproduction of heavy 
quarks at NNLO where the heavy quark mass dependence 
is fully retained. Our results form part of the NNLO de- 
scription of heavy quark pair production relevant for the 
NNLO analysis of ongoing experiments at the Tevatron 
and planned experiments at the LHC. 

Our analytical results arc presented in factorized 
forms. For the divergent parts, the squared matrix el- 
ements are given in terms of the Laurent series expan- 
sion of the corresponding LO and NLO contributions ex- 
panded up to 0(e) and 0(e 2 ), respectively. In this way 
we could transfer parts of the finite part of the squared 
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amplitudes to the coefficient functions of the pole terms. 
After this we found that the remaining parts of the finite 
contribution could be further factorized, partly in terms 
of the corresponding LO and NLO pieces, and, for the 
box graphs, in terms of factorizing forms as described in 
Sec. [V] Writing our analytical results in factorized forms 
led to a reduction of the length of the original output by 
a factor of 400. To the best of our knowledge these nice 
factorization properties of the squared amplitude were 
not known before. It would be interesting to find out the 
underlying reason for this factorization. 

The present paper deals with unpolarized quarks in the 
initial and final states. Since our results for the matrix 
elements in [l3[ contain the full spin information of the 
process, an extension to the polarized case with polariza- 
tion in the initial state and/or in the final state including 
spin correlations is not difficult. 

The present calculation constitutes a first step in ob- 
taining the full NNLO corrections to the heavy quark pro- 
duction processes in QCD. A further next step is to pro- 
vide one-loop squared results for gluon-initiated heavy 
quark pair production. Work on the gluon-initiated 
channel is in progress. 



Analytical results in electronic format for all the terms 
in Eq. (|4.1[) . including the (t <-> u) symmetric terms ex- 
plicitly written out, as well as combined full results, are 
readily available [23| . 
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APPENDIX 

Here we present the expressions for the terms Q m appearing in Equation (|5. 10[) . 



Qi = 

Q 2 
Q 3 



Q 5 

Qe 
Qi 



8m 2 (s/t - z 2 /{s0 2 )) - B{ 0) 2(m A s + D)/t- B^'2z u /0 A + C^ 1 (2D + s 2 + 2tz t + 2m 2 z 2 /0 2 ) 



(0), 



C { ° ] 2t(s + 2T) + C y i'\2D + s 2 + 2t 2 ) - D y 2 >t(2D 



f(0) 



(o), 

2 



2t A 



ID 



2/t - B[ a) jt + c[ 0) Tz 2 /D - C^htT/D + Cf\l + t 2 /D) - D { 2 0) t{l + t 2 /D) , 

4(2m 2 z t - D)/(sl3 2 t) + B[ 0) 2T/t + Bf ] 2z t / {s/3 2 ) + c[ 0) (z t /f3 2 + (m 2 st - t 3 )/D) 
+ Ci 0) 2t 3 /D + C { ° ] st 2 /D - D ( 2 0) S t 3 /D , 

8m 2 z u /{sl3 2 ) + Bfhm 2 + Bf ) 2m 2 z 2 /{sl3 2 ) + cf ] {2tT + m 2 z 2 //3 2 ) 
- Cf ] 2m 2 t + Cf ) (m 2 s + 2t 2 ) - D { 2 0) t(m 2 S + 2t 2 )] /D , 

16m 2 z u /(s/3 2 ) + B { °hm 2 + B ( ° ] Am 2 z 2 / (s/3 2 ) + c[ 0) 2(2tT + m 2 z 2 /f3 2 ) 



C { °hm 2 t + C^>2(m 2 s + 2t 2 ) - D 2 u >2t(m 2 s + 2t 2 ) /D , 



(0), 



16m 2 /(s/3 2 ) + B[ { '>2 - B^2/P 2 - C[ (4:t(D + m 2 t) + s 2 T//3 2 + Am 2 t 2 /(3 2 )/D 
+ cf ] 2tTz 2 /D - cf ] z 2 {D + t 2 )/D + Df ] tz 2 (D + t 2 )/D , 

8m 2 (s/t - 4 - 5z t /(sf3 2 )) - B { " ] 2(2D/t - 3m 2 + u) + Bf ] 2{m 2 s + 6m 2 t - su)/(s(3 2 ) 

+ C[ 0) (2m 2 s + WtT+(m 2 + s)z 2 /(3 2 ) ~ cfht^m 2 + z 2 ) + cf\s{hm 2 + z 2 ) + 10t 2 ) 
- D { 2 0) t(s(5m 2 + z 2 ) + m 2 )] /D , 



(A.l) 
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Qs = 

Q 9 = 
Qio — 



Ql2 
Ql3 

Qu 

QlB 
Ql6 



Ql7 
Ql8 



8m 2 (D/t - tz 2 /{s(3 2 )) - B[ 0) 2T(2D/t - s) + B { 5 0) 2(D + tz t /(3 2 ) 

- C[ 0) s{m 2 s - t 2 - tzt - tz t /p 2 - t 2 (m 2 s - t 2 )/D) - C^ ] 2stT{l + st/D) 

+ Cf ) s(m 2 s + 2t 2 + st 3 /D) - D^st^s + 2t 2 + st 3 /D) , 

-4(T/i + z t /(s(3 2 )) + B[ a) 2T/t + Bfhzt/isp 2 ) + c{ 0) (z t /(3 2 + t(m 2 s - t 2 )/D) + cf } 2t 3 /D 
+ Cf ) st 2 /D-D$ ) st 3 /D , 

8m 2 D - BfhD Zt + B<g ] 2tD - C[ 0) st(m 2 s + Am 2 t + t 2 ) + cf } 2i 2 (m 2 s - t 2 ) 



+ c{ 0) st(m 2 s -t 2 )- D$"st 2 (m 2 s - t 2 ) /D , 



(0) o+ 2, 



On = \%Dz u /{s[3 2 ) + B[ a) 2D + B ( £ >) 2z 2 D/{sf3 2 ) - C[ 0) s{m 2 s - t 2 - z 2 D/{s(i 2 )) - C^hst 2 



-Ci 0) s 2 t + D' v 's 2 t 



(O) D 2 .2 



8m 2 z t /(s/3 2 ) + S| 0) 2T - Bf ] 2{D - tz t )/(sf3 2 ) + c{ 0) sz t ((2m 2 - s)/(s(3 2 ) + t 2 / D) 

- ci 0) 2st 2 T/D - cflsHT/D + D ( 2 0) s 2 t 2 T/D , 

Cfszt + Cf 2t 2 + cf ) st-D { 2 a) st 2 , 
8m 2 {s/u - z 2u /(s{3 2 )) - B{°j2(m 2 s + D)/u - B^hzt/p 2 + {2D + s 2 + 2uz u + 2m 2 z 2u /p 2 ) 

- C£>2u(s + 2U) + ci 0) (2D + s 2 + 2u 2 ) - D { ^u(2D + s 2 + 2u 2 )] /D , 

2/u - B^/u + c[ 0) Uz 2u /D - cilhuU/D + cf\l + u 2 /D) - £><°>u(l + u 2 /D) , 
~8m 2 (s/u - z u /(s(3 2 )) - B™2(D - m 2 t)/u - ^ 0) 2(i - m 2 z 2u /( S (3 2 )) 

+ Ci 0) {s 2 - tz u + D//3 2 - 2m 2 uz 2u /(sf3 2 )) - Cfj2u{m 2 + z 2u ) 
+ Cf ] (D + s 2 - tu) - D ( °Ju(D + s 2 - tu)] /D , 

16m 2 /(s/3 2 ) + B$2 - Bfh/p 2 - c[ 0) {4u(D + m 2 u) + s 2 U//3 2 + Am 2 u 2 //3 2 )/D + C^2uUz 2u /D 

- cf ] z 2u {D + u 2 )/D + D^uz 2u (D + u 2 )/D , 

8m 2 (4s/u - 1 - 5z u /(sf3 2 )) - B[ j2(5D - 3m 2 t + tu)/u - B^2(4t - 5m 2 z 2u /(s[3 2 )) 

+ c[ 0) (As 2 + 2(4s + 5u)U + 5m 2 z 2u /l3 2 ) - C^2u(5m 2 + Az 2u ) + ci 0) {5m 2 s + Asz 2u + 10u 2 ) 

- D ( ^u{5m 2 s + 4sz 2u + 10u 2 )] /D . 
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